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This contribution gives an overview on primordial non-Gaussianities from a theoretical 
perspective. After presenting a general formalism to describe nonlinear cosmological per- 
turbations, several classes of models, illustrated with examples, are discussed: multi-field 
inflation with non-standard Lagrangians, modulaton fields, curvaton fields. In the latter 
case, a special emphasis is put on the isocurvature perturbations, which could leave a spe- 
cific signature in non-Gaussianities. 

§1. Introduction 

A potentially promising probe of the early Universe, which has been studied very 
actively in the last few years, is the non-Gaussianity of the primordial perturbations. 
Whereas the simplest models of inflation, based on a single field with standard kinetic 
term produce undetectable levels of non-Gaussianity, a significant amount of non- 
Gaussianity can be produced in scenarios with i) non-standard kinetic terms; ii) 
multiple fields; iii) a non standard vacuum; iv) a breakdown of slow-roll evolution. 

In this contribution, after summarizing the formalism to describe primordial 
non-Gaussianities, we review three categories of models that could produce de- 
tectable non-Gaussianities. First, we discuss inflationary models with generalized 
Lagrangians involving multiple scalar fields and non standard kinetic terms. Sec- 
ond, we consider non-Gaussianity generated by modulaton fields, or light scalar 
fields which are spectator fields during inflation but can affect some cosmological 
transition. Finally, we present some aspects of non-Gaussianity arising from curva- 
ton models, in particular the issue of isocurvature non-Gaussianity. 

§2. Nonlinear cosmological perturbations 

In this section, we first present a geometrical description of the nonlinear cos- 
mological perturbations. We then summarize the traditional statistical description 
of primordial perturbations, which are used to relate early Universe models with 
cosmological observations. 

2.1. Covariant approach 

Instead of the traditional metric-based approach, 1 ) our discussion below is based 
on a more geometrical approach. 2 ) Let us consider a spacetime with metric g a b and 
some perfect fluid characterized by its energy density p, its pressure P and its four- 
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velocity u a . The corresponding energy momentum-tensor is given by 

T ab = pu a u b + P(g ab + u a u b ). (2-1) 
Let us also introduce the expansion along the fluid worldlines, 

= V a u a , (2-2) 



and the integrated expansion 



N = l I dr0 , (2-3) 



where r is the proper time defined along the fluid worldlines. By noting that 0/3 
corresponds to the Hubble parameter H in a homogeneous and isotropic spacetime, 
one can interpret 0/3, in the general case, as a local Hubble parameter and a\ oc = e N 
as a local scale factor, N being the local number of e-folds. 

The conservation law for the energy-momentum tensor, X7 a T a b = 0, implies 3 )' 4 ) 
that the covector 

satisfies the relation 

(oS£ „ ( „ = -_!_(v o P^V„p) , (2-5) 

where a dot denotes a Lie derivative along u a , which is equivalent to an ordinary 
derivative for scalar quantities (e.g. p = u a V a p). 

If w = P/p is constant, the above covector is a total gradient and can be written 

as 

Ca = V a N+ 1 hip . (2-6) 
3(1 + if) 

On scales larger than the Hubble radius, the above definition is equivalent to the 
non-linear curvature perturbation on uniform density hypersurfaces 5 * 1 

( = SN- ( P H d -l = 5N + l - [ P -^-. (2-7) 

The above equation is simply the integrated version of (2-4). 

The covector Q a , or the corresponding scalar quantity Q, can be defined for 
the global cosmological fluid or for any of the individual cosmological fluids (this 
approach can also be extended to the case of interacting fluids 6 )). 

2.2. Power spectrum and higher order correlation functions 

We now discuss the statistical properties of the cosmological perturbation Q. We 
first define the power spectrum 

(C kl Ck 2 ) = (2vr) 3 <5( 3 )(k 1 + k 2 )P c (A;i) , (2-8) 
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where the Fourier modes are defined by 

C k = (2vr) 3 J d 3 x e" ik - x C(x) . (2-9) 

Beyond the Gaussian case, the first quantity of interest is the three-point func- 
tion, or its Fourier transform, called the bispectrum and defined by 

(C kl Ck 2 Ck 3 ) = (2n) 3 5^C£ k ^ B d k ^ k ^ k ^- ( 2 - 10 ) 

i 

Equivalently, one often uses the so-called /nl parameter, which can be defined in 
general by 

B c (h, k 2 , k 3 ) = ^hUh, k 2 , k 3 ) [P f (fc!)P c (fe) + P c (k 2 )P c (k 3 ) + P f (fc 3 )P c (fci)] . 

(2-H) 

Similarly, the Fourier transform of the connected four-point function defines the 
trispectrum, according to 

(Ck 1 Ck 2 Ck 3 Ck 4 ) c ^(2vr) 3 <5( 3 )(^k l )T c (k 1 ,k 2 ,k3,k4). (2-12) 

i 

The various correlation functions defined above can in principle be measured, or at 
least constrained, by observations of CMB fluctuations or large scale structure. 

In many models, the final perturbation £ depends on the fluctuations of one or 
several scalar fields generated during inflation and the so-called 5N- formalism 7 )' 8 ) 
provides a powerful method to evaluate, at least formally, the primordial non- 
Gaussianity generated on large scales. 9 ) The underlying idea is to describe, on scales 
larger than the Hubble radius, the non-linear evolution of perturbations generated 
during inflation in terms of the perturbed expansion from an initial hypersurface 
(usually taken at Hubble crossing during inflation) up to a final uniform-density hy- 
persurface (usually during the radiation-dominated era) . Using the Taylor expansion 
of the number of e- folds given as a function of the initial values of the scalar fields, 

one finds, 9 )' 10 ) in Fourier space, 

(c kl Ck 2 Ck 3 ) = Y, n j n ,j n M^IMM 3 ) + 

UK 

l - N I N J N KL {6^i 1 6ipi 2 (6ip K *6<p L ) k3 ) + perms. 

IJKL 

(2T4) 

From this expression, one can distinguish two types of non-Gaussianity, depending 
on whether the first line or the second line dominates. 
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In the first case, non-Gaussianities arise from the three-point function of the 
scalar field(s), as models with non-standard kinetic terms. 11 )~ 13 ) This leads to a 
specific shape of non-Gaussianity, called equilateral because the dominant contri- 
bution comes from configurations where the three wavevectors have similar length 
k\ ~ hi ~ &3. 

In the second case, non-Gaussianities arise from the nonlinear dependence of N 
of the scalar field(s). Assuming quasi-Gaussian scalar field fluctations, one finds that 
(2-14) leads to a bispectrum of the form (2-11) with 

6 NjNjN" 

s /nl " {n k n«y ■ (245) 

This corresponds to another shape of non-Gaussianity, usually called local or squeezed, 
for which the dominant contribution comes from configurations where the three 
wavevectors form a squeezed triangle. 

The present observational constraints 14 ) on these two main types of non-Gaussianity 

are 

-10 < 4)° cal) < 74 (95% CL), -214 < /^ uil) < 266 (95% CL) . (246) 

Note that other types on non-Gaussianity can be produced, such as the folded shape 
(which peaks at h% ~ hi + hi) arising from a non-standard vacuum. 13 ) 

Extending the Taylor expansion (2-13) up to third order, one can compute in a 
similar way the trispectrum. 15 ) For local non-Gaussianity, the trispectrum can be 
written in the form 16 ) 

T c (ki,k 2 ,k3,k4) = t NL [P(h 13 )P(h 3 )P(h 4 ) + 11 perms] 
54 

+^5nl [P(.k 2 )P(h 3 )P(h 4 ) + 3 perms] , (2-17) 

with 

N LJ N IK N J N K 2hN IJK N I N J N K 
TNL = (N L Ntf ' 9nl = Va {NlNL? (248) 

and where h\ 3 = |ki + ka|. The present constraints on these parameters, assuming 
the data do not contain any isocurvature contribution, are 17 ) 

-7.4 < W- 5 g NL < 8.2 (95% CL), -0.6 < 10" 4 rjv L < 3.3 (95% CL) . 

The next section will be devoted to models with several inflatons described by 
a very general Lagrangian, with multi-field DBI inflation as a concrete illustration. 
Later, we will consider in turn two types of scenarios leading to local type non- 
Gaussianity: modulatons and curvatons. 

§3. Generalized multi-field inflation 



Although it is obviously simpler to deal with a single inflaton, high energy physics 
models usually predict the existence of many scalar fields. Multi-field inflation should 
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thus be considered seriously. Moreover, many models involve kinetic terms that are 
not standard. It is therefore useful to develop a formalism that is as general as 
possible and try to tackle multi-field models described by an action of the form 



S = I d A X\f^g 



R +P(X IJ ^ K ) 



16vrG 



(3-1) 



where P is an arbitrary function of N scalar fields and their associated N(N + l)/2 
kinetic terms 

X IJ = -Iv^vV- (3-2) 

This includes multi-field models of the form 

P(X,4> I ) = X-V(</> 1 ), X = G U X IJ , (3-3) 

where Gjj is an arbitrary metric in field space, as well as k- inflation models, 18 ) which 
depend on very general single- field Lagrangians. 

The expansion up to second order in the linear perturbations of the action (3-1) is 
useful to obtain the classical equations of motion for the perturbations and to calcu- 
late the spectra of the primordial perturbations generated during inflation. Working 
for convenience with the scalar field perturbations Q 1 defined in the spatially flat 
gauge, the second order action can be written in the compact form 19 ) 

S(2) = \J dtd 3 xa 3 [(P<jj> + 2P <MJ>t<IK> X MK ) Q T Q J 

-PKij^diQfdjQ 3 -M KL Q K Q L + 2Q ki Q k Q i '\ (3-4) 

where the mass matrix Mkl an d the mixing matrix is Qki can be expressed in 
terms of the X IJ , of P and its derivatives. 

An interesting example, which combines non-standard kinetic terms with a pos- 
sibly multi-dimensional inflaton space, is multi-field DBI (Dirac-Born-Infeld) infla- 
tion. In these models, motivated by string theory, inflation is due to the motion of 
a L>3-brane in an internal six-dimensional compact space. 20 )' 21 ) The dynamics of 
the brane is governed by the Dirac-Born-Infeld Lagrangian, hence the name of these 
models. 

After some appropriate reparametrizations, the Lagrangian can be written in 
the form 

(^/det^ + fGud^d^J) -lj- Vtf 1 ), (3-5) 



where the potential arises from the brane's interactions with bulk fields or other 
branes. This Lagrangian can be written explicitly in the form (3-1) upon using 

det(5£ + / GudPtfdvt 3 ) = 1 - 2fG u X IJ + 4f 2 X [ /xj ] 

-8f 3 X [ fxjX$ + Wf'xfxjxZx? (3-6) 

where the field indices are lowered by the field metric Gjj, which corresponds to the 
metric of the internal compact space, and the brackets denote antisymmetrization 
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over the indices. Note that the dilaton and the various form fields are ignored in the 
Lagrangian (3-5), but they can also be included in the analysis of the cosmological 
perturbations generated by these models. 22 ) 

As in other multi-field inflationary models, it can be convenient 23 ) to decom- 
pose the perturbations into a so-called (instantaneous) adiabatic mode, along the 
background velocity in field space, and (instantaneous) entropic modes, which are 
orthogonal to the adiabatic direction. Focussing, for simplicity, on the two-field case, 
where there is a single entropic degree of freedom, one can decompose the scalar field 
perturbations as 

Q 1 = Q<r e a + Qsel , (3-7) 

where the adiabatic vector e£ and entropic vector el are normalized (via the field 
space metric Gu). The perturbations generated during inflation can then be deter- 
mined by using the standard techniques, which gives 24 ) 

H 2 H 2 
V Q°* - 4^2 > V Qs* - ( 3 - 8 ) 

(the subscript * here indicates that the corresponding quantity is evaluated at sound 
horizon crossing kc s = aH) . For small c s , the entropic modes are thus amplified with 
respect to the adiabatic modes. 

Since we are in a multi-field scenario, the curvature perturbation can evolve after 
sound horizon crossing, and the spectrum of the final curvature perturbation, which 
is probed by cosmological observations, can be formally written as 

Vn = {l + Tl s )Vn, = (1 + ^s)^ = (1 + ( 3 ' 9 ) 

where T-jis quantifies the transfer from the entropic into the adiabatic modes. 

Let us now discuss non-Gaussianities in multi-field DBI inflation. The three- 
point correlation functions of the scalar fields can be computed from the third order 
action, which is given, in the small sound speed limit, by 19 )' 24 ) 

5 (3) = J dtd 3 X [(Qa) 3 +C 2 MQs) 2 

~2^& [Q°( dQ,r f ~ ^QviPQ,) 2 + 2c 2 Q s dQ a dQ s )\ I . (3-10) 

The contribution from the scalar field three-point functions to the coefficient /nl is 
found to be given by 

f(3) _ 35_ 

i08 C 2 (i + r 2 5 )' 

which is similar to the single- field DBI result, but with a suppression due to the 
transfer between the entropic and adiabatic modes. 

Interestingly, multi-field DBI inflation could also produce a local non-Gaussianity 
in addition to the equilateral one. 25 ) Finally, let us mention that the trispectrum in 
multi-field DBI inflation has also been computed. 26 ) 
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§4. Modulatons 

Significant non-Gaussianity can arise when a cosmological transition in the his- 
tory of the Universe depends on some light scalar field, which has previoulsy acquired 
some fluctuations during the inflationary phase. Consequently, in different regions of 
the Universe where the value of the scalar field is slightly different, the cosmological 
transition and the subsequent cosmological evolution will differ. In this way, the 
fluctuations of the scalar field, which we will call a modulaton, are converted into 
curvature fluctuations. 

4.1. Modulated reheating 

A typical example is the modulated reheating scenario 27 )' 28 * 1 where the decay rate 
of the inflaton, r, depends on a modulaton a. 

A simple way to compute the curvature perturbation is to calculate the number 
of e-folds between some initial time ij during inflation, when the scale of interest 
crossed out the Hubble radius, and some final time tf. For simplicity, let us assume 
that, just after the end of inflation at time t e , the inflaton behaves like pressureless 
matter (as is the case for a quadratic potential) until it decays instantaneously at 
the time t d characterized by H d = r. At the decay, the energy density is thus 
Pd = p e exp[— 3(N d — N e )] and is transferred into radiation, so that, at a subsequent 
time tf, one gets 

p f = pte-Wf-"*) = Pe e-Wf- N <)-( N f- N *). (4-1) 

Using the relation r = H d = H f exp[2(N f - N d )} to eliminate (N f - N d ) in (4-1), 
we finally obtain 

N f = iV e --m^---ln — . (4-2) 
3 p e 6 Hf 

This implies that the nonlinear curvature perturbation can be simply expressed as 

< = < M -iln(q?>), (4-3) 

where Cinf represents the contribution to the curvature perturbation from the inflaton 
fluctuations. At linear level, this leads to the curvature power spectrum 

n - + ^ (£) 2 n„, - p c „ + JL (£) 2 (f ) 2 . (44, 

By expanding (4-3) up to second and third orders in 5a*, one can easily determine the 
bispectrum and trispectrum for the curvature perturbation. Using the expressions 
(2-15) and (2-18), one finds that the associated nonlinear parameters are given by 



rr»\ w2 

and 



/ NL = 5(l- T ^)H 2 , (4-5) 



TNL = !/£ L H- 5NL = f (2-3^ + ^)^ (1-0) 
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where = 1 — (T-'( inf /'P() represents the fraction of the curvature power spectrum 
due to the modulaton. 

4.2. Modulated trapping 

One can also envisage the possibility that a modulaton field affects the cosmolog- 
ical evolution during inflation. This is the case in the modulated trapping scenario, 29 ^ 
which relies on the resonant production of particles during inflation. 30 ) In this model, 
the inflaton <j> is coupled to other fields, for example to some fermions ip via with the 
interaction Lagrangian £; nt = Xcptpip. 

If, during inflaton, the effective mass of ip, m e ff = m — X(f> becomes zero, this 
triggers a burst of production of these particles. Even if these particles are quickly 
diluted by the exansion, their backreation will affect the evolution of the inflaton, 
governed by the equation of motion 

4> + 3 H <p + V' (0) = N\{^) = An* (— \ G{t-U). (4-7) 

Indeed, the term on the right hand side induces a temporary slow-down of the 
inflaton, which leads to a slightly longer phase of inflation. This brief trapping of 
the inflaton thus manifests itself as an increment of the number of e-folds until the 
end of inflation: 

N = N std {0) + Z\iV trapping (4-8) 

Let now assume that this trapping depends on some modulaton, for example via 
the coupling between the inflaton and the particles, and occurs well after the fluctua- 
tions of the modulaton (on observable cosmological scales) have been generated. This 
is in contrast with other scenarios 30 ) -33 * 1 where the trapping occurs approximately 
when cosmological scales exit the Hubble radius, which leads to special features in 
the CMB spectrum as well as specific non-Gaussianity. 34 - 1 

Using the Taylor expansion of (4-8), where only the second term AN depends 
on the modulaton a: 

C = SN= + • • • + AN^Sa + \AN^ba 2 + \AN, ac M* (4-9) 

(higher order derivatives with respect to the inflaton are ignored, because they 
give negligible non-Gaussianities), one can compute the power spectrum and non- 
Gaussianity of the curvature perturbation, generated by the modulated trapping 
scenario. 29 ^ According to (2-15), the corresponding non-linearity parameter for the 
bispectrum is given by 

Similarly, the nonlinear coefficients of the trispectrum are 

(AN aa ) 2 =3 36 2 25 AN „„„ ^ 3 

= (AN„r " _ ^ /NL ' 9NL - 54 (ZiV^ " ( } 
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The most interesting situation occurs when the coupling A depends directly on 
the modulaton, which leads to the nonlinear parameters 



3 + 2 



AA^\ „• 
A' 2 



5NL 



2e 2 /3 2 



1 AA" A 2 A'" 
1 + 6— r + 4- 



A' 2 



3A' 3 



=■3 



(4-12) 



where /3 = Max(zA</>)/|(^>*| cannot exceed 1. 

If A depends only linearly on a, the first expression reduces to /nl = 3£" 2 /(2e/3), 
which shows that it is quite easy to obtain a detectable level of non-Gaussianity in 
this scenario: for example, one gets f^L — 55 with (3 = 0.01 and E = 1. Moreover, 
there is a specific relation between tnl and (/nl which could be confronted with obser- 
vations if these quantities can be measured, and thus distinguish this scenario from 
other scenarios leading to different relations between the nonlinear coefficients. 35 ) 

If the future cosmological data point to the existence of a significant amount 
of local non-Gaussianity, the modulated trapping scenario would thus represent a 
viable model, together with the modulated reheating or the curvaton scenario which 
we examine in the next section. 



§5. Curvatons and isocurvature perturbations 

The last example that we consider in this contribution is the curvaton scenario, 36 ) 
or more precisely in the mixed curvaton and inflaton version 37 ) where the inflaton 
fluctuations are also taken into account. The curvaton is a weakly coupled scalar 
field, a, which is light relative to the Hubble rate during inflation, and hence acquires 
Gaussian fluctuations with an almost scale-invariant spectrum. After inflation the 
Hubble rate drops and eventually the curvaton becomes non-relativistic so that its 
energy density grows with respect to that of radiation, until it decays. 

Many aspects of the curvaton scenario have been studied in the literature. Here, 
we wish to focus our attention on isocurvature perturbations that can be generated in 
this type of scenario, and their non-Gaussianity. Isocurvature non-Gaussianity, which 
has been investigated recently in several works, 38 ) -44 ) could indeed be distinguished 
from the usual adiabatic non-Gaussianity and thus open a new window on the early 
Universe, if ever detected. 

As a preliminary step we present a general formalism that computes systemati- 
cally the evolution of the nonlinear perturbations of various fluids through a decay 
transition. We then apply this formalism to a scenario with a curvaton fluid, ra- 
diation and cold dark matter CDM and compute the adiabatic and isocurvature 
perturbations, up to third order. 

5.1. Evolution of the perturbations due to the decay of some species 

We now consider a very general setting where several cosmological fluids coexist, 
each of them characterized by the nonlinear curvature perturbation 

a = W + -^ rln^, 

3(1 + 71?^) PA 



(5T) 
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as follows from the definition (2-6). We wish to compute the curvature perturbations 
after the decay of one these fluids, denoted a, which will later correspond to the 
curvaton. 

In the sudden decay approximation, which we adopt here, the decay takes place 
on the hypersurface characterized by = r a , where is the Hubble parameter 
at the decay and r a is the decay rate of a. Since H depends only on the total energy 
density, the decay hypersurface is a hypersurface of uniform total energy density, 
with 5Nd = £, where £ is the global curvature perturbation. The equality between 
the sum of all energy densities, before the decay and after the decay, thus reads 

£p A _ e 3(i+^)(a-C) = ^ = £p fl+e 3(i4- fl )(C B+ -C> > ( 5 .2) 

A B 

where the subscripts — and + denote quantities defined, respectively, before and after 
the transition. In the above formula, we have used the non-linear energy densities of 
the individual fluids, which can be expressed in terms of their curvature perturbation 
(a by inverting the expression (5-1). 

Expanding the first equality in (5-2) up to third order, one finds 



A 

with the coefficients 



Ca- + % (Ca- - C) 2 + $ (Ca- - cf 



(5-3) 



p A = 3(l + w A ), A A = 2p ft A = (l + w A )ft A , n = Y,&Ai (5-4) 

n A 

where the abundance parameters are defined just before the decay: ft A = p A -/ p&ec&y 
Note that, although the global perturbation ( appears on both sides of (5-3), this 
relation can be used iteratively in order to determine, order by order, the expression 
of C in terms of all the Ca~, up to third order. 

Just after the decay, the energy density of the fluid a is transferred into the one 
or several of the remaining fluids. Introducing the relative branching ratios j Aa , this 
means that the energy density for any species A, just after the decay of a, is simply 
given by 

PA+ = PA- + lAaPa- (5-5) 

This relation, which is fully non-linear, can be reexpressed, upon using (5T), in the 
form 

e M( A+ -0 = (j _ / A ) e / ? A(CA--C) + f Ae M<v--0 j ( 5 . 6 ) 
where the parameter 

f A = "lAattg ^_ 7 j 
ft A + 1 Act ft a 

represents the fraction of the fluid A that has been created by the decay. 
Expanding (5-6) up to third order, and using (5-3), one gets 



C4+ = £t 



B 
A 



Cb- + f (Cb- - C) 2 + § (Cb- - C) 3 
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-^(Ca+-C) 2 -^(Ca+-C) s 



Pa 
2 



6 



with the coefficients 



T A A = f A [l ~ 
Tf = f A U- 
T A C = f A (l- 



Pa 
Ax 
Pa 

Pa_ 

Pa 



Xa + (1 - fA) 
K + fA lT A 

Ac, C^A,a. 



(5-8) 

(5-9) 
(5-10) 
(5-11) 



Finally, after using (5-3) again, one finds that (5-8) yields the full expression, up to 
third order, of all the post-decay curvature perturbations in terms of the pre-decay 
curvature perturbations: 43 ^' 44 ) 



C4+ = E T As- + E u% c (b-Cc~ + E Va cd (b~Cc-Cd-, 

B B,C B,C,D 



(5-12) 



with 



tjBC - 1 
A = o 



E Pe t a {Seb - Xb)(S ec - A c ) - Pa(Tab - Xb)(T ac - X c ) 



and 



Vf CD = -\Y,PeT A e^>eb - Xb)XfPf(Sfc - X c ){Sfd - X D ) 



E,F 



+ g E Pe t ae{5eb - Xb){5ec - Xc)(Sed - Ad) 
e 



U% D - \ E /3sA s (fe c - X c )(5ed - A D ) 



~Pa(Tab - Xb) 

-qP A (T A b - Xb){T ac - Xc)(T A d - Ad) • 

The above expression thus provides a systematic computation of the post-decay 
curvature perturbations for all fluids in a very general setting. For scenarios with 
several decay transitions, the perturbations can be obtained by combining the various 
expressions of the type (5-12) for each transition. 

5.2. Mixed curvaton and inflaton scenario 

We now apply the general formalism presented above to a scenario involving a 
curvaton a, behaving as a pressureless fluid, in addition to radiation (r) and CDM 
(c), which can lead to isocurvature perturbations. 45 ^' 46 ^ The formula (5-12) allows 
us to compute, in terms of the pre-decay perturbations, the perturbations Q r and Q c 
after the decay, or equivalently the adiabatic perturbation, which coincides with C, r 
deep in the radiation era, and the CDM isocurvature perturbation 



Sc — 3(Cc — Cr)- 



(5-13) 
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For simplicity, we restrict our analysis to the situation where 

Cc- = Cr- = Cinf , (544) 

by assuming that the CDM and radiation perturbations, before the curvaton decay, 
depend only on the inflaton fluctuations. 

The curvaton fluid isocurvature perturbation before the decay, S a , can be easily 
related to the curvaton field fluctuations in the case of a quadradic potential. Indeed, 
writing the (non-linear) energy density of the oscillating curvaton defined on the 
spatially flat hyper surf aces, characterized by SN = Q r when the curvaton is still 
subdominant: 

p a = m 2 a 2 = m 2 (a + 5a) 2 = p a e 3 ^ a ~ Cr) = p a e StT , (5-15) 
leads to the relation 



5a \ 



2 



e s " = [1 + yj . (5-16) 

Expanding this expression up to third order, and using the conservation of 5a /a in 
a quadratic potential, we obtain 

Sa = S-±S 2 + ±S 3 , (5-17) 

where the quantity 

S = 2— (5-18) 

is Gaussian. 

Using the general expressions (5-12), one finds that the primordial curvature 
perturbation is given by 

Cr = Cinf + ZlS + l ~Z 2 S 2 + ^ 3 ^ 3 , (5- 19) 



with 



f r f 4r r 2 \ 

Zx = V Z2 = Is v 3 " 8r + T"¥j' (5 ' 20) 

r 2 SGr 3 24r 2 4r 2 48r Ihr 18 \ , 

where the parameter 

C - : n 7rg ]Q (5-22) 

1 - (1 - 7r CT )J/ ( T 

can be interpreted as the efficiency of the energy transfer from the curvaton into 
radiation (£ = 1 if the curvaton decays only into radiation, i.e. j ra = 1), and 



Primordial non-Gaussianities 



13 



For the primordial isocurvature perturbation, one obtains 



S c = Sl S + \s 2 S 2 + ^s 3 S 3 , (5-24) 



with 



1 / „ , . 2r 3 4r 



2 



si = f c -r, s 2 = ^ ( 3/ c (l - 2/ c ) + ^ - + 8r 2 - 3rJ , (5-25) 

l,2/o r 2 /6r 3 24r 2 4r 2 48r 15r 18 \ 

«3 = - g/c (3 " 4/c) - — + - - — - + 64r + — - 36 j . (5-26) 

From these expressions, one can determine the power spectrum, the bispectrum and 
the trispectrum, which can be probed by observations. 

The power spectrum for the total curvature perturbation follows from the linear 
part and is given by 

^=^ + J^ = S " 1 J^. ( 5 - 27 ) 

where E represents the fraction of the power spectrum due to the curvaton con- 
tribution. The power spectrum for the isocurvature fluctuations is, according to 
Eq. (5-24), 

n -IT—,- ( 

Both curvature and isocurvature perturbations depend on the curvaton fluctuations 
and are therefore correlated, with the correlation coefficient: 



■P Sc = (f c -r) 2 V s . (5-28) 



e f y/S, £/ = sgn(/ c - r). (5-29) 



In the pure curvaton limit (E ~ 1), adiabatic and isocurvature perturbations are 
either fully correlated, if ef > 0, or fully anti-correlated, if Ef < 0. In the opposite 
limit (E <g 1), the correlation vanishes. For intermediate values of E, the correlation 
is only partial, as can also be obtained in multifield inflation. 47 ) 
The isocurvature-to-adiabatic ratio 

a = ^ £ = 9 ( 1 "7) 2 H ' t 5 ' 30 ) 

is strongly constrained by cosmological observations, the precise limits depending on 
the assumed level of correlation between the isocurvature and adiabatic perturbations 
(since the impact of isocurvature perturbations on the observable power spectrum 
depends crucially on this correlation 48 ) ) . In terms of the parameter a = a/(l + a), 
the limits (based on WMAP+BAO+SN data) given in 14) are 

a < 0.064 (95%CL), ai < 0.0037 (95%CL) , (5-31) 

respectively for the uncorrelated case (E = 0) and for the fully correlated case 
(E = 1). According to (5-30), the observational constraint a <C 1 can be satisfied 
if \ f c — r\ <C r (which includes the case f c = 1 with r ~ 1) or if E <C 1, i.e. the 
curvaton contribution to the observed power spectrum is very small. 
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5.3. Adiabatic and isocurvature non-Gaussianities 

Since we now deal with two observable quantities, namely adiabatic and isocur- 
vature perturbations, the definition of the bispectrum can be extended to include 
both types of perturbations. In our particular case, where there is only one degree 
of freedom, S, at the nonlinear level, one can show that the generalized bispectra 
(with indices / = {£, S}) are of the form 43 ) 

B IJK (h,k 2 ,h) = 6^P s (fc2)P s (fc 3 ) + 6jf P^P^ + b^P^P^k,) 

(5-32) 

with 

^r J L K = ^(i)^ )t (5-33) 

where ArL = z 2 , N^. = s 2 , AT/L = z±, = s\, respectively. 

Recalling that the usual, purely adiabatic, /jvl is proportional to the bispectrum 
of ( divided by the square of the power spectrum, one defines the analogs of /tvl by 
dividing the coefficients b 1 ^ 3 ^ by the square of the ratio Pf/Pg = z 2 /E, i.e. 

;i,jk _ 6 f i,JK _ j JK 

f NL =^NL =^f b NL ■ ( 5 - 34 ) 

Taking into account the fact that the last two indices can be permuted, this leads to 
six different coefficients, explicitly given by the expressions 

fCCC _ Z 2 „ 2 K,(S _ SlZ 2 2 r5,CC _ S 2 j-,2 /k qrrA 

J JVL — 9 " ' JiVL — ^3 " ' J NL — 9 " ' 
z l z l z l 

?<,SS _ 5 1 Z 2 „ 2 ?S,CS _ g l s 2 „ 2 ?5,55 _ s l$2 „ 2 /r o fi \ 

Jatl — 4 — , JjVL — 3 - , JjVL — 4 — • ^O-ODJ 
z 1 z 1 z 1 

The same analysis applies to the trispectra that combine adiabatic and isocur- 
vature perturbations, leading to the generalized parameters 44 ) 



U,KL = 4)^(1)^) .3 -I.JKL = Ur JKL = N^N^ 
NL — z 6 " ' &NL — 25 gNL ~ z 6 



(5-37) 

where = z% and = S3. Taking into account the symmetries under per- 
mutations of the indices, one finds, for two observables (I = {(,S}), 9 different 
parameters t^^ L and 8 parameters g N J L KL '■ 

An interesting question is whether one can find significant non-Gaussianities, 
while satisfying the bound on the isocurvature spectrum. As mentioned earlier, 
this isocurvature constraint can be satisfied with E ~ 1 if f c and r are sufficiently 
close. In this case, one finds that the purely adiabatic non-Gaussianity dominates. 
But in the alternative situation where E <C 1, one finds that, with respect to the 
purely adiabatic non-Gaussianity, the purely and mixed isocurvature ones are either 
enhanced by constant factors, if f c <C r <C 1, or much more strongly enhanced with 
powers of (3/ c /r), if r <C f c <C 1. 
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§6. Conclusions 

As this contribution has tried to illustrate with a few explicit examples, the 
detection of primordial non-Gaussianities would have dramatic consequences. 

First, since the simplest inflationary models, based on a single field in slow roll, 
predict a negligible amount of primordial non-Gaussianities, these models would 
have to be replaced with more elaborate models, involving several scalar fields, non 
standard kinetic terms or other features. 

Second, since the measurement of non-Gaussianities contains potentially a lot 
of information, in particular concerning their shape, one could hope to discriminate 
between different categories of models, which would otherwise appear degenerate in 
their predictions of the power spectrum. 

Even the simplest shape of non-Gaussianity, the local shape, can hide sur- 
prisingly rich variations if perturbations are generated by several scalar fields and 
if isocurvature perturbations survive. In such situation, purely adiabatic, purely 
isocurvature and mixed non-Gaussianities could coexist and the hierarchy between 
their amplitudes would provide invaluable information. Since isocurvature fluctua- 
tions are usually associated with the generation of dark matter and baryon asym- 
metry in the Universe, non-Gaussianity from isocurvature fluctuations, if detected 
in the future, would give us a lot of insight into the nature of dark matter, the 
mechanism of baryogenesis, and therefore into high energy physics. 
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